Pion Compton scattering and bremsstrahlung 

Goran Fald J] and Ulla Tengblacfl 

Department of nuclear and particle physics, 
Uppsala University, Box 535, S-751 21 Uppsala, Sweden 

Abstract 

The pion-polarizability functions are structure functions of pion-Compton scattering. They can 
be assessed in high-energy pion-nucleus bremsstrahlung reactions, 7r~ + A — > tt~ + 7 + A. We 
present numerical expectations for pion-nucleus bremsstrahlung cross sections in the Coulomb 
region, i. e. the small-angle region where the nuclear scattering is dominated by the Coulomb inter- 
action. We investigate the prospects of measuring the polarizability functions for pion-Compton 
cm. energies from threshold up to 1 GeV. A meson-exchange model is used for the pion-Compton 
amplitude. 
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I. INTRODUCTION 



The cross-section distribution for high-energy pionic bremsstrahlung 
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is at small-momentum transfers to the nucleus dominated by the one-photon-exchange con- 
tribution. Hence, by measuring cross-section distributions for pion-nucleus bremsstrahlung 
at small momentum transfers one gains at the same time information about the cross-section 
distribution for pion-Compton scattering 



7 + 7T 



7 + 7r~ 



the photon of the initial state being the virtual photon exchanged between pion and nucleus. 
In particular, it becomes possible to determine pion polarizabilities. 

The theoretical underpinnings of the pion-bremsstrahlung reaction were investigated by 
Gal'perin et al. and an experiment was susequently performed by Antipov et al. 2], which 
showed that pion polarizabilities could indeed be determined. The COMPASS experiment 
at CERN J3( is a refinement of the Antipov experiment. 

In ref. |4|] a mathematical expression for the bremsstrahlung cross section was derived 
assuming a pion-Compton amplitude that in additon to the Born contributions contained 
also contributions from the a, p, and ai, exchanges. The kinematics of the nuclear reaction 
was defined through 

7T-(pi) + A(p) -> 7T-(p 2 ) + T(?2) + A(p% (1) 

and the kinematics of the related pion-Compton reaction through 



(Pi) + 7(?l) -»• T (P2)+7(?2), 



(2) 



with qi = p — p' . In terms of these kinematic variables the cross-section distribution can be 
written as {4] 
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This expression is valid when the transverse momenta of the emerging pion and photon 
are much smaller than their longitudinal momenta. The parameter a is the fine-structure 
constant, m n the pion mass, x the ratio 

PI V ^ 

= qij_ ■ q2_i_, and q m in the minimum-longitudinal- momentum transfer 

_ m|_ x 

qmm ~2E 1 'l-x [) 
At high energies q m in is extremely small. 

The functions A(x, q|j_) and -B(x, q|j_) characterize the pion-Compton amplitude. In the 
Born approximation 

A(x,ql ± ) = l, (6) 
B(x,ql L ) = 0; (7) 

and generally 



A{x, qL) = 1 - T7T~~~ T f ^tt^ ) M*, <&-). ( 8 ) 

qL) = 7uT~ \ ( ql± t X 2 X l A 2 (X, qL), (9) 
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with Ai(x, q|_|_) and \2{x, qfj_) generalized-polarizability functions. Their mathematical ex- 
pressions, in the meson-exchange model, are given in 4J. 

The structure of the cross-section distribution is mainly determined by the inelastic- 
Coulomb factor 

n 2 

9l± (10) 
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which vanishes when the transverse-momentum transfer q l± to the nucleus vanishes. This is 
the classic Primakoff factor, which exhibits a maximum at q±± = q m in- When the momentum 
transfer q%± becomes much larger than q m in the Coulomb contribution to the bremsstrahlung 
cross section becomes small and the hadronic contribution dominates [5] . This contribution 
is not considered in the present paper. 



The Mandelstam kinematic variables of the pion-Compton scattering can in our ap- 
plication be approximated as follows; 



1 



s — mz 



x(l — x) 



2,22 

q 2± + x m n 



1 — x 



2,22 

q 2± + x m w 



(11) 



u — mz 



2,22 

q 2± + x m w 



4 



II. KINEMATIC AL CONSIDERATIONS 



We are interested in studying the pion-Compton-scattering process. In order to do so in 
a transparent way we need to make specific cuts in the pion-nucleus-bremsstrahlung phase 
space. To this end a few kinematical observations might be helpful. 

Let us first introduce, instead of the pion-Compton cm. energy and the photon-transverse 
momentum, the corresponding dimensionless variables 

w = \fs/m n (12) 
V = |q2±IA?V, (13) 



with w > 1 and y > 0. The first entry of Eq. (flTl) . defining the pion-Compton cm. energy 
squared, can then be reformulated as 

y 2 = x[(l-x)w 2 - 1]. (14) 

Since y 2 is positive we conclude that for a given fixed value of w the variable x must be in 
the interval 

< x < x max = 1 . (15) 

w 2 

The maximal value of y occurs at the midpoint of this interval, 

w 1 

2/max= -(1-^)- (16) 

We may also rewrite Eq. (ri~4l as 

y 2 = x(x max - x)w 2 , (17) 

showing directly that the photon-transverse momentum vanishes at the kinematic end-points 
x = and x = x max . 

For a fixed value of x, on the other hand, it follows from the positivity of y 2 of Eq. ( [T4l 
that the pion-Compton cm. energy is bounded from below 

w 2 >l/(l-x). (18) 

The boundary value obtains in the forward direction, where the photon-transverse momen- 
tum vanishes. 
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In Fig.l the photon-transverse momentum y(x, w) is plotted for a few values of the pion- 
Compton cm. energy w = s/s/m n . As an illustration consider ^Js = 2m.,,, a mass at which 
the threshold approximation to the polarizability functions should be acceptable. In this case 
the maximal value of a; is 0.75 and the maximal value of the photon-transverse momentum 
q 2 ± is 3m w /A. 
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FIG. 1: Photon-transverse momentum y = q2±/m 7V as a function of x for fixed values of the 
Compton mass w = ^/s/m w = 2, 3, 4, 5. The threshold value, w = 1, corresponds to the point 
y = x = 0. 
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III. LOW-MASS-COMPTON REGION 



Hard-pion bremsstrahlung offers a unique possibility of determining pion-threshold po- 
larizabilities. However, this will succeed only if we restrain the cm. mass in the pion-photon 
system to masses well below the p-meson mass, so as to avoid the s-channel p-exchange 
contribution to the Compton amplitude. 

Hence, we first concentrate our efforts on the region s -C m 2 , where we can confidently 
approximate the structure functions A(x, qfjj and B(x, qfe ± ) of Eqs.© and (J9j) by their 
near-threshold values, with 

\ x {x, q 2 ± ) = — * K + fa) = 0.0006, (19) 
a 

X 2 (x,4 ± ) = fa = -0.0131, (20) 
a 

and a n and the standard pion-electric and pion-magnetic polarizabilities. These values 
given are obtained in the meson-exchange model when fixing the cx-meson coupling so that 
the chiral Lagrangian value |g of A2 is reproduced. Our present knowledge of the pion 
polarizabilities is reviewed in |7|. 

The cross-section distribution (j3J) contains angular dependent terms, through \i = qi± ■ 
q2j_. Measured variables are however q 2 j_ and p2_i_, but since, p 2 j_ = qu_ — q 2 j_, with 
Qi± q2±, the angular variations of q 2 j_ and qu_ are independent. Thus, the integrations 
over the angles of q 2 ^ and qi^ in Eq.Q are unrestricted and we may replace /i 2 by its 
angular average, which is a half, and at the same make the replacement 

dV±d 2 g 2 ± = 7i 2 dq 2 1± dq 2 2± . (21) 

After integration over angles the cross-section distribution factorizes into factors that 
depend either on qi± and x, or on q 2 ± and x. The integration over the inelastic- Coulomb 
factor (ITU]) is an integration over qi±, but there is an x dependence in q m i n , 

qmm ~2E 1 1-x [ ] 

It is tempting to integrate over all possible momentum transfers qi±, but that is not possible 
since the integral diverges. The reason is that we have neglected pion and nuclear form 
factors. Also, for large values of qi± the hadronic-bremsstrahlung contribution dominates 
the Coulomb contribution. Consequently, we must make a cut so as to avoid the hadronic 



contribution. The best approach is to make a cut in the same way for all x values. This is 
achieved by cutting off the integral at r^q min on the low-momentum side and at r + q min on 
the high-momentum side of the Coulomb peak. From fllOp we derive 

Pf = P(r-,r + ) = / dgl ± q ^ (23) 

This integral is obviously independent of x. The precise values of r_ and r + to be chosen 
can only be decided after inspection of the experimental distribution. 

Next we integrate the cross-section distribution Eq.Q over transverse-photon momenta 
in the domain < q\ x < q\^ max - Since the polarizability functions Ai and A2 entering 
the structure functions A(x,c$±) and B(x,qfe ± ) of Eqs.(|HD an d & are now constants the 
integration is straightforward, and we get 

dcr(x,q 2 ± m ax) _ 4:Z 2 a 3 1- x 
dx m 2 x 

with 



P f F(z), (24) 



z(x,q 2±max ) = q ^ (25) 

and the distribution function 

2r 2 / I 1 \ 

F{z) = u{z) + [ ri v(z)X 1 + -w(z)X,j . (26) 

The x-dependent prefactor in Eq. (|24p is simply the ratio of the final state momenta, 
(1 — x)/x — E2/U2, and the singularity at x ~ is the well-known soft-photon-radiation 
singularity. 

The functions u(z), v(z), and w(z), of Eq. fT26l) are all elementary and given in the Ap- 
pendix. They are graphed in Fig. 2a. For small values of the argument z, or transverse 
momenta such that q%± ma x ^ a;2m ^ we have 

u{z) w v(z) w w(z) w z. (27) 

For large values of the argument z, or transverse momenta such that q\^ max 3> x 2 ml, we 
have instead 

2 

m(z) w -, v(z) 2, iw(z) ~ 1. (28) 

In particular we note for large values of z the ratio u(z)/w(z) = 2/3. 

An alternative to integrating over transverse-photon momenta is to integrate over Comp- 
ton masses y/s in the domain m 2 , < s < s max . This is simply achieved by making the 
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FIG. 2: Plots of u, v, and w as a) functions of z; and b) functions of x for yjs max = 3.75™,^. In b) 
we have plotted f/10 instead of v. 



integration limit q2±max depend on s max . According to Eq. ffTTj) . or (TBI) , the relation be- 
tween q2±max and s max allows us to write 

Z\X) S max ) n 7r~ Tri^max 2<) > 



2 



with 

%max = 1 • (30) 

°max 

At the same time as being the definition of the variable z, Eq. (129|) displays the functional 
relation between q2± ma x and s maa; . 

We now specialize to y/s max = 3.75m n which is the cut for the threshold region employed 
by the COMPASS collaboration. This implies, via Eq. (!30|) . a maximal x- value of x max = 
0.93. In Fig. 2b the three distribution functions u, v/10, and w are plotted as functions of x 
for this value of s max . From this figure we draw the following conclusions. 

Near the upper end of the x-distribution the three functions u, v, and w, may be ap- 
proximated by their threshold values, Eq. fl27|) . This is a consequence of Eq. (1291) . which 
demonstrates that even though the factor s max /m^ = 14 is large, when x is sufficiently 
close to x max the variable z will be small. This observation leads to a simplification of the 
cross-section-distribution function F(z) of Eq. (j26]) . 



F(z) 



1 + ^(^-1^; 

m± 



(31) 



Numerically, we expect A 2 — §Ai ~ —0.0134 and small, but enters multiplied by a large factor, 
Smax/^nii and therefore makes an important contribution to the cross-section distribution. 



In the mid-range x-region the value of z is large, due to the large value of s max /rn%, and 
we may here approximate the functions v and w by their asymptotic expressions, Eq. (T28l) . 
This allows us to simplify the expression for the function F(z), yielding the approximate 
formula 



F\x^ s max t 



U\X) s max ) 



3s r 



3a; 



(x\ 2 + 5A1 



(32) 



Ami ""' ' 2(1 ■ ./')"■ 2 
The function u(z), representing the Born approximation, must be evaluated exactly as we 
are looking for small deviations from it. The numerical value of the coefficient multiplying 
the middle Ai term in Eq.( I3"2"]) is large. Nevertheless, this term has essentially no effect on 
the cross section since according to Eq. lTiTIj) Ai is proportional to the sum a n + j3 n which 
vanishes, or is very small compared with f3 n . The important term is the A2 term. Over most 
of the x-range its strength is fixed by Eq. fl32l . But this approximation cannot be used near 
the upper end-point since it would there overestimate the polarizability contribution by a 
factor of 3/2, as a comparison with expression (J3Tj) would demonstrate. 
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FIG. 3: Plots of a) the distribution function F(z(x, s max )) of Eq. ()26p in the Born approximation 
(dashed line), and in the full calculation (solid line); and b) the ratio of the two distribution 
functions. 

In Fig. 3a the cross-section-distribution function F(z) of Eq. (|26p is graphed in the full cal- 
culation and in the Born approximation, i.e., when polarizability parameters Ai = A2 = 0. 
In Fig. 3b the ratio of the two distributions is plotted. Clearly, the polarizabilities are impor- 
tant only in a region near the end-point of the x- distibution, and therefore the experimental 
efforts are concentrated on this region Q]. 
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In the COMPASS experiment the apparatus is blind in a small region around the forward 
direction. The transverse momenta of the final state pions and photons must be larger than 
Q±cut to be detected. As a result the distribution function F(z) of Eq. (12oT) must be replaced 
by 

F(z) - F(z) - F(z ) (33) 



where from Eq.fl25 



z (x) 



x 2 m1 



(34) 



In the COMPASS experiment q± cu t = 15 MeV/c and Zq is small except when x is near 
x = 0. The subtraction procedure of Eq.( l33l works unless there are momentum transfers 
<?2± < q±cut for which s > s max . This will not happen as long as x remains in the domain 



Q±cut 



< x < 1 



m 



Q±cut 



mi 



mi 



(35) 



which is the case in the COMPASS application. The effect of the cut on the cross section is 
illustrated in Fig. 4. 




FIG. 4: The cross-section-distribution function when small-momentum transfers are cut out ac- 
cording to Eq. (|33p . The dashed line is the distribution function without the cut, F(z), and the 
solid line with the cut, F{z) — F{zq). 
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IV. RHO-MASS-COMPTON REGION 



At large Compton masses the polarizability functions are complex and depend in a compli- 
cated manner on the varables s and x. As before, after angular integration, the cross-section 
distribution factorizes into one factor dependening on q 11 _ and x, and another one depending 
on s and x. 

The integration over the Coulomb-peak factor is unchanged and yields the factor Pf of 
Eq. (l23l) . As the Compton energy s increases and we move closer to the end-point of the 
x distribution, the minimum-momentum transfer, Eq.([5]), increases. The maximal value is 
obtained at the end-point itself, which is located at x = 1 — m^/s, and where q m in — s/2E\. 
At a Compton mass equals the rho mass and an energy E\ = 190 GeV this means q m in = 1.5 
MeV/c. Away from x ~ 1 the value is much smaller. 

The angular integration of the right hand side of Eq. ([3]) is also unchanged. However, it 
is convenient to replace the variable q\± by s, variables which are related by Eq. ffTTj) . The 
cross-section distribution obtained is 

dcr 4Z 2 a 3 „ . 

i^ = j^w P ' GM - (36) 

The distribution function G(x, s) is given by the expresssion 

G(x, s) = |(1 - Y(x, s))A(x, s)\ 2 + \Y(x, s)A(x, s) + B(x, s)\ 2 , (37) 
with s related to via Eq.tfTTl) and 

q_2j_ + x z m^ 1 — x s — m4 

with the obvious restriction < Y < 1. 

It is instructive to look at the structure functions as functions of x and s. From Eqs (jSJ) 
and (Ej) we derive the expressions 

A(x,s) = l-l^£(t=^j X 1 (x,s), (39) 

B(x,s) = \ A 2 (x,,). (40) 

In the Born approximation A = 1 and B = 0, as stated in Eqs. ([6]) and (JJJ). The Born 
contribution, which is by far the dominant contribution, is in the cross-section distribution 
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(|37|) multiplied by the polynomial factor 

p(Y) = 1-2Y + 2Y 2 , (41) 
with Y = Y(x, s). Since < Y < 1 it follows that 

~ < POO < 1. (42) 

For a fixed value of s, the maximal value of the pre-factor, p(Y) = 1, is obtained at the 
boundary points x = and x = 1 — rn^/s. The minimum value, p(F) = 1/2, is obtained 
at a; = 1 — 2m^/(s + m^). Thus, for large values of s/m 2 there is rapid variation in the 
kinematical factor p(Y) near the upper end-point of the x region. In the remaining part of 
phase-space p(Y) is smoothly varying and near to one, since there the variable Y(x, s) itself 
will be quite small. This is all illustrated in Fig. 5. 
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FIG. 5: A contour plot of the kinematic factor p(Y(x, s)), or equivalently, the Born approximation 
to the distribution function G(x, s). Compton mass s in units of m 2 . The dashed line represents the 
kinematic boundary, where p(Y) = 1; and the dash-dotted line the minimum, where p(Y) = 1/2. 

The pion structure functions A and B depend on the Mandelstam variables s, t, and u, of 
the pion-Compton scattering, or equivalently x and s. In the region we are considering the 
Born terms dominate and the polarizability terms mainly come in through their interference 
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with the Born terms. In Fig. 6a we plot the quantity \A — 1| 2 . From Eq. (1391) it is clear that 
this is a measure of the strength of the Ai term. The graph shows that the Ai term is small, 
except for x- values near x = 0, that is grows with s and exhibits a structure related to the 
p- and ai-meson exchanges. 




FIG. 6: Three-dimensional plots of the functions a) \A — 1| 2 and b) B 2 . Compton mass s in units 
of m 2 . 

In Fig. 6b we plot the quantity \B\ 2 which measures the strength of the A2 term. The 
graph demonstrates that A2 is important in a region of phase space, essentially disjoint 
from that of Ai, namely the region near the upper boundary line of x- values. Its strength 
grows with increasing value of s. The structure related to the p- and ai-meson exchanges is 
visible but less pronounced. It is A2 that contains the cr-meson exchange term. According 
to Eq. (1371) the interference between B and the Born term has an extra factor Y(x, s). This 
factor results in a sharper structure in the interference term than seen in B alone. 

Fig. 7 illustrates the behaviour of the dynamic factor G(x, s) of the integrand (13T1) . The 
plotted quantity is the ratio G(x,s)/Gb(x,s), i.e. the ratio of G(x,s) to its Born approx- 
imation Gb{x,s) = p{Y). We notice that in the overwhelming part of phase space the 
Born approximation is reasonably accurate. The exceptions are at small x-values where 
the polarizability function Ai causes substantial deviations from unity, and at x-values near 
the upper-boundary line where the polarizability function A2 contributes to the structure. 
Unfortunately, in the latter region precise measurements are complicated by the sharp struc- 
ture in the distribution function G(s,x). 

In Fig. 8a the distribution function G(x, s) is graphed as a function of s for fixed x = 0.2. 
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FIG. 7: The ratio of the cross-section-distribution function G(s,x) to its Born approximation 
Gb(x,s) = p(Y). The Compton mass s is measured in units of m^, and the ratio x = ufy/Ei. 

The dependence on the pion-polarizability functions is strong, with a first dip at the p-meson 
mass. As pointed out above this structure is essentially caused by the A x contribution. This 
is seen already in Eqs. (1391) and fHOj) where the Ai contribution is multiplied by the factor 
(1 — x) and the A2 contribution by the factor x. In Fig. 8b the distribution function G(x, s) 
is instead graphed as a function of x for fixed s = 30 m^. This Compton mass is right on 
top of the p-meson mass. Again we observe that the deviation from the Born contribution is 
larger for small x-values than for large x-values. Near the upper end-point there is a narrow 
dip, a structure mainly caused be the ^-dependence. 
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FIG. 8: The cross-section-distribution function G(s,x) as a) a function of s, in units of m 
x = 0.2; and as b) a function of x for s = 30 m^. 
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V. SUMMARY 



Hard bremsstrahlung in high-energy pion-nucleus collisions in the Coulomb region is 
induced by single photon excitation. Therefore, this reaction can be used to extract infor- 
mation about pion-Compton scattering. 

The pion-Compton amplitude is dominated by Born terms that give the exact amplitude 
for point-like pions. The non-pointlike structure is described by two independent polariz- 
ability terms. We describe them in a meson-exchange model with a, p, and ai, mesons. The 
polarizability contributions are small but are important to study. At threshold in the pion- 
Compton system Ai oc (a n + f3 n ) is predicted to vanish and A 2 oc f3 n predicted to be small. 
We investigate, through the pion-bremsstrahlung mechanism, the polarizability functions 
for pion-Compton cm. energies s below 1 GeV. 

In that part of phase space where most of the final-state energy is carried by the photon 
(x ~ 1), the contribution from Ai is negligible. The contribution from X 2 is small but 
its strength increases with increasing s. Some structure caused by the p-and ai-meson 
exchanges can be seen. At the pion-Compton threshold the contribution from the a-exchange 
term dominates. Experimental efforts to measure the polarizability contributions are at the 
moment concentrated to this region. 

In that part of phase space where most of the final-state energy is carried by the pion 
(x m 1), the contribution from X 2 is negligible. The contribution from Ai is clearly seen, and 
exhibits a strong variation due to the p-and ai-meson exchanges. There is no contribution 
from a exchange in A x . 

Our analysis is valid at high energies when the transverse momenta are small compared 
with the longitudinal momenta. In addition the momentum transfer to the nucleus must be 
in the Coulomb region, i.e. very small. 



17 



VI. ACKNOWLEDGEMENTS 



We would like to thank Jan Friedrich for valuable discussions and for information about 
the COMPASS experiment, and Bengt Karlsson for help with the presentation. 



18 



VII. APPENDIX 



In Sect. 3 we integrate the cross-section distribution Eq.Q over the photon-transverse 
momenta q 2 ±, from the lower limit q 2 ± = up to the upper limit g 2 _i_ = l2_Lmax- Then, we 
encounter the following elementary integrals; for the integration of the Born contribution 

u(z) - z [ q ^ max dq ^ ( x2m * Y [i _ 2x2m l c i± \ 

+ ^—r,i (43) 



(1 + ^)2 3(1 + 2)3 

for the integration of the Ai contribution 

viz) - z r^ max ll- 2x2m ^ <1 ^ ) 

= ^±^-ln(l + ,); (44) 

and for the integration of the A2 contribution 

[Z 2 2± ma x dq 2 ± ( x 2 ml \ f q 2 ± 1 

w \ z ) = z / ~2 — — , — 2 — 2 M 1 ~ 2 2 i — T~ ( 



1 + z 

The parameter z is defined as 



(45) 
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